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In the pseudogap regime of one-dimensional incommensu- 
rate Peierls systems, fluctuations of the phase of the order 
parameter prohibit the emergence of long-range order and 
generate a finite correlation length. For classical phase fluc- 
tuations, we present exact results for the average electronic 
density of states, the mean localization length, the electronic 
specific heat and the spin susceptibility at low temperatures. 
Our results for the susceptibility give a good fit to experimen- 
tal data. 

PACS numbers: 71.23.-k, 02.50.Ey, 71.10.Pm 



Continuous symmetries in one-dimensional electronic 
systems are not spontaneously broken at any finite tem- 
perature T. The mean-field prediction of a finite critical 
temperature T^ F is incorrect in this case. The experi- 
mentally observed Peierls transition at a finite temper- 
ature Tj? D <C in many quasi one-dimensional con- 
ductors is due to weak interchain-coupling which triggers 
a crossover to three-dimensionality. In the intermediate 
temperature regime < T < T^ F , the physical prop- 
erties of Peierls chains are dominated by one-dimensional 
order parameter fluctuations Q. This is the so-called 
pseudogap regime where mean-field theory is not even 
qualitatively correct. 

In this work, we shall present an exact solution of an 
effective model for the low temperature thermodynamics 
of incommensurate Peierls chains. For incommensurate 
chains, the order parameter A(x) is complex, so that in 
the pseudogap regime the generalized Ginzburg-Landau 
potential has the form of a "Mexican hat" [0] . It is then 
a good approximation to ignore amplitude fluctuations 
of A(x) = |A(x)|e rf ^- ) and focus on the gapless fluctu- 
ations of the phase #(x). At long wavelengths and low 
energies, the thermodynamics of phase fluctuations can 
be described by a classical Hamiltonian which is formally 
identical with the kinetic energy of a superfluid pj , 



H,<> = — m*n, 
2 



dxv 2 s {x) , v s (x) 



2m* 



(1) 



Here, m* is the effective mass of the electrons, L is the 
length of the chain, and the one-dimensional density n s 
measures the stiffness of the system with respect to long- 
wavelength distortions of the phase of the order param- 
eter gj. A two-dimensional analog of Eq. ([]]) has been 
used by Emery and Kivelson || to explain the pseudo- 
gap behavior of underdoped high-temperature supercon- 
ductors. In one dimension, the problem is much sim- 
pler, because there is no Kosterlitz-Thouless transition 



and the thermodynamics of the phase variable is triv- 
ial. However, the calculation of the electronic properties 
amounts to solving a one-dimensional random problem 
with colored noise. Usually, problems of this type can- 
not be solved exactly At low energies, the electronic 
degrees of freedom can be described by the Hamiltonian 
of the so-called fluctuating gap model (FGM) 



H c \ = —ivpd x a3 + A(x)a + + A*(x)a- 



(2) 



where vp is the Fermi velocity, cr^ are the usual Pauli ma- 
trices, with o± = h(ci ± i&2)- In the pseudogap regime, 
the gapped amplitude fluctuations are frozen out, so that 
we may set A(x) — A s e l ^ x \ where A s is determined by 
the local minimum of the generalized Ginzburg-Landau 
functional. From Eq. ([!]), it is then easy to show that 



(A(x)A*(x')) 



-\x-x'\/£ 



(3) 



where (...) denotes the thermodynamic average with re- 
spect to the Hamiltonian H$ given in Eq. (fil), and the 
order parameter correlation length is £ = n s /(2m*T). In 
this work, we shall calculate the average electronic den- 
sity of states (DOS) p(u) = (Tr£(o; - H cX )) of the model 
defined via Eqs. ([l])-^) exactly for arbitrary £. 

Previously, the DOS of the FGM has been calculated 
assuming a Gaussian distribution of A (a;) with covariancc 
given by Eq. @ Although in this case the prob- 

lem is not exactly solvable ]8||| , a sophisticated algorithm 
has been developed JjJ which produces an expression for 
p{uj) which is reasonably close to the exact numerical 
result for Gaussian disorder |^]. However, as explained 
above, the assumption of a Gaussian distribution of A(x) 
centered at A = is rather unphysical in the pseudogap 
regime. It is therefore not surprising that in this regime 
the true behavior of p(u>) (to be discussed below) is very 
different from the DOS for Gaussian disorder. 

The electronic contribution to the thermodynamic 
properties of our system can be obtained from the 
disorder-averaged free energy 



-T 



dx{p(x,u))Yn(l + e~ w / T ) . (4) 



The local DOS p(x, u>) can be expressed as p(x,u>) — 
— 7r _1 ImTrt/(a;, x, u) + i0 + ), where the Green function 
G(x, x' , u) + i0 + ) satisfies 



[cj — H c \] G(x, x' , uj + i0 + ) = S(x — x')(Jq 



(5) 



Here, <jq is the 2x2 unit matrix. For periodic boundary 
conditions, the average {p{x,uj)) is independent of a; and 
can be identified with the average DOS p(u>). 
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Following Ref. , we eliminate the phase of the order 
parameter A(x) = A^e'*^ via a gauge transformation, 



(6) 



The transformed Green function Q satisfies an equation 
of the form (0) , but with H e \ replaced by 



H c \ = -iv F d x (T 3 + V(x)a + A s ai 



(J) 



where V(x) = ^f-dx^x). Eq. (^J) is a chiral transfor- 
mation that eliminates the phase of A(x) in favor of 
a forward scattering random potential V{x). The lo- 
cal DOS is invariant under this transformation, so that 
we may replace Q — > Q in all expressions involving the 
DOS. The crucial point is now that with H§ given by Eq. 
(|l|), the probability distribution of V{x) is determined by 
Gaussian white noise, with zero average and covariance 
(F(x)F(x')) = v%{2g)- 1 8{x - x'). Due to the Gaussian 
white noise statistics of V(x), the average DOS of our 
model can be calculated exactly in the thermodynamic 
limit. Several methods of obtaining the exact p(u>) are 
available. Actually, Eq. (]?]) is a special case of the class of 
random Hamiltonians discussed by Hayn and Mertsching 
jnj , who calculated the average DOS by means of the su- 
persymmetry method jig] . Alternatively, the DOS can 
be calculated within the phase formalism ||. In Ref. 
jl3| a modification of this formalism [|||] is used to di- 
rectly obtain the integrated Green function T(lo) defined 
via 9 w r(a;) = Tr (G(x, x, uu + i0 + )) from the solution of a 
Fokker-Planck equation. For L — > oo, only the stationary 
solution is needed, and we obtain 

r(w) = £-\uj) - inAf(u) = npoAJ^M/I-M , (8) 

where po = (nvp)^ 1 is the DOS for A = 0, and Ii V {g) is 
a modified Bessel function with imaginary index iv. The 
dimensionless parameters g and v are 



4A s g 
v F 



2n s A s 
m*VF T 



4u£ 
vf 
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A, 



(9) 



Note, that the imaginary part of is proportional 

to the integrated average DOS JV(oj) which satisfies 
d u J\f(uj) = p(uj) while, according to Thouless |Q, the 
real part of r(w) can be identified with the inverse mean 
localization length £~ 1 (uj), i.e. the Lyapunov exponent 
||. Using a Wronski relation for I-i v (g) (ll| we get 



PoVf sinh(7r^) 



4 < \h„{g)\ 2 ' 

For the inverse mean localization length we get 

£-V) = ^|-ln|/ w ( 3 )| . 
v F og 



(10) 



(11) 



We now discuss the behavior of the average DOS. Be- 
cause p{oj) is an even function of w, we restrict ourselves 



to us > 0. Using 15 |J i;/ (0)| 2 = ( 7riy )^ 1 sinh(7w) one 
easily verifies that Af(w) ~ p^oj for g — > 0, so that in 
this limit we recover the result for free electrons with lin- 
earized energy dispersion. While for small g, the leading 
corrections can be calculated perturbatively in powers of 
g, in the pseudogap regime j > 1, the behavior of the 
average DOS is quite complicated. It is convenient to 
measure frequencies in units of A s and to express Eq. 
( pp| ) in terms of the Bessel function Ji V (ig) with imag- 
inary index and argument, using Ii U (g) = e" 71 '/ 2 Ji V (ig) 
fll5|. Defining Q — uj/A s — v/g, we may write 



Po 



\J igiI1 {ig)X 



(12) 



In Fig. U, we show a graph of Eq. ([l2|) for several values of 
g. For a more quantitative analysis, we use the uniform 
asymptotic expansion of Ji g a{ig) for large g and fixed Co 
JlB]] which reveals three different regimes: First of all, for 
1 — Co > <? -2 / 3 (i.e. for frequencies sufficiently far below 
A s ), the average DOS in the pseudogap regime g ^> 1 
can be approximated by 



p(w)/po « 2g(l - u 2 ) 1 ' 2 exp[-2 fl Q(<D)] 



-27rgu)~\ 



arccos(u;) 



(13) 



where Q(u>) = (1 — w 2 ) 1 / 2 — w arccos(w) . In particular, 
for small Q, we may expand Q(u>) ~ 1 



T 0J 



1 ,-,2 



or, so 



that 



p(u))/p a w 2it ge~ 29 cosh[-K gu) e' 9ul , gQ 6 < 1 . (14) 

Hence, for uo — 0, the DOS is exponentially small, 
p(0)/po ~ 2ixge~ 29 . As shown in Fig. || such a strong 
suppression of the DOS at the Fermi energy is a unique 
feature for classical phase fluctuations, which is neither 
reproduced within the Born approximation H (which 
predicts p(0) cx ff _1 ), nor for Gaussian disorder [fl (where 
p(0) cx <?~ p , with /i w 0.64). The approximation (H3) 
breaks down when 1 — a) becomes comparable with 

Note that Q(l - e) - =nr e3/2 for e < 1, so that 
gQ{u>) = 0(1) when Eq. ( |l3| ) ceases to be valid. In this 
case, we have to go back to our exact result (112) which 



implies for 



11 < g- 2/3 « l 



(15) 



Here, a x = 2- 4 / 3 7r- 1 c 2 /c? w 0.7306 and a 2 = 
2 2/3 [3(c 2 /ci) 2 - ci/c 2 ] « 0.3534, with ci = Ai(0) = 
[3 2 / 3 r(2/3)]- 1 and c 2 = -Ai'(O) = [3 1 / 3 r( 1/3)]" 1 , 
where Ai(x) is the Airy function. From Eq. (|l5|), we 
conclude that, to leading order in g ^> 1, the average 
DOS exhibits a maximum precisely at w = A s , with a 
hight that diverges as g 1/3 cx £ 1/3 cx T" 1 / 3 for T -> 0. 
Finally, for a) — 1 > g~ 2 ^ 3 , our exact result ( |l2|) reduces 
to the well known expression for the DOS in the presence 
of a static gap, p(oj)/po w u/\fu} 2 — 1. At cD — 1 w g~ 2 ^ 3 , 
this expression smoothly matches with the parabola jl^). 
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In Fig. ||, we show the exact inverse localization length 
£ _1 (u;) given in Eq. ( [ill ) f° r several values of g. For 
g > 1 we obtain the following approximations: I (ui) w 
(A s /V)(l-u> 2 ) 1/2 for l-o) > <T 2 / 3 ; for |w-l| < .$T 2/3 
we find ~ {^s/v F )[a 3 g~ 1/3 - (Sg)" 1 ], with a 3 = 

2- 2 / 3 c 2 /ci re 0.4592; finally, forw-1 > 5~ 2 / 3 the leading 
behavior is « (A s /w f )[2 3 (w 2 - l)]" 1 . 

Let us now consider the electronic contribution to the 
free energy F c \ defined in Eq. (Q). For the FGM with a 
linearized energy dispersion, Eq. (Q) is ultraviolet diver- 
gent, because then the DOS approaches a constant for 
\uj\ — ► oo. However, physical quantities involve deriva- 
tives of F e i, which at low temperatures depend only 
on the low-energy part of the spectrum and are finite. 
For convenience, we regularize Eq. (Q) by subtracting 
from F c \ the free energy F^°° for an infinite correla- 
tion length, where the gap is static. After an integration 
by parts, we express the integral in Eq. (^) in terms of a 
fermionic Matsubara sum and obtain 



F r , - Fz 



2LA S T 



E 

n=0 



, (16) 



where Q n — n(2n + 1)T/A S . For large Q n , the term 
in the square bracket vanishes as cl>,7 2 , so that the sum 
converges. In the pseudogap regime g ^> 1, we may use 
the uniform asymptotic expansion of Iga n (g) for large 
g fl5|| to obtain an expansion of Eq. ( jig ) in powers of 
g^ 1 cx cx T. For T <C A s the leading terms are 



F r ] - f: 



^— oo 



L 

16? 



12< 



+ o(r 2 ) 



(17) 



The physical interpretation of this result is simple: Be- 
cause £ is roughly the size of domains where the order 
parameter is spatially constant, the prefactor L/£ in Eq. 
( |l7| ) can be identified with the number of locally ordered 
domains in a system of size L. At distances of the order 
of £, the phase fluctuations distort the order parameter, 
which leads to an increase of the energy. In the limit 
A s £/t>_F — * oo, the energy scale associated with a twist 
in the order parameter is set by A s . For finite £ this en- 
ergy scale decreases, because the time £,/vf it takes for 
electrons to propagate over the distance £ is finite. This 
gives rise to the second term in Eq. (|l7j). We emphasize 
that our exact result ( |l6| ) gives the change in the free 
energy due to phase fluctuations for arbitrary £. 

The low-temperature behavior of the specific heat 
C c i = —Td 2 F c \/dT 2 can be calculated analytically. 
Keeping in mind that £ = n s /(2m*T), we see that the 
leading contribution to C c \ is due to the first correction 
term (involving the energy vf/0 in Eq. (|17|), 



C el ~ (n 2 /24){n /n s ) 2 p LT , 



(18) 



where no = rn'vpj'K, and we have used the fact that 
the contribution from F^°° is exponentially small due 
to the static gap. Thus, in the pseudogap regime, the 



electronic specific heat of Peierls chains is linear in T, 
just as the specific heat for non-interacting electrons in 



one dimension, C, 



(0) 



-,<o) 



\PqLT. Note that C cX /C y cl 
^(no/n s ) 2 for T — > 0. In general, we expect that n s /no 
is a number of the order of unity for T <C T c MF [fil|, so 
that C c i/C^ = 0(1). In the same regime, we find from 



Eq. (|14|) that p(0)/ Po ~ 4^exp[-f i. e. the 

average DOS at the Fermi energy is exponentially small 
(see the dashed line in Fig. |J). 

Given p(u>), we may also calculate the spin susceptibil- 
ity x = T- 1 J °° dup{w) cosh" 2 (w/2T). A graph of X as 
a function of T/T C MF is shown in Fig. | (solid line). The 
low-temperature behavior can again be calculated ana- 
lytically. If g 3> 1 but n s /no < |, we find \ ~ l/o(0) (as- 
suming now s = 2 for spin degeneracy 0] ) . On the other 
hand, for n s /no > j the frequency integral is dominated 
by a new saddle point at u> = cosr, 
Using Eq. (|l3|), we obtain 



where r = 

8 n s 



x/xo ~ 2(2 7 r) 1 /V(A I ,/T) 3 / 2 exp[-A I ./T] 



(19) 



where A r = ^^As and xo = 2po is the susceptibility 
of free electrons. n s /no > i implies r < f , so that at 
low temperatures the ratio X /2p(0) is exponentially large, 
x/2p(0) cx exp[(l - sinr)A s /(rT)]. Our graph of x(T) 
in Fig. ^ is quite similar to the corresponding graph 
given by Lee, Rice, and Anderson [||. Note, however, 
that these authors assumed a real order parameter and 
an exponentially large correlation length at low temper- 
atures. Because incommensurate Peierls chains are char- 
acterized by a complex order parameter and a correlation 
length that diverges only as a power law, £ cx T _1 , the 
agreement between the theoryof Ref . ||] and experiments 
for incommensurate chains |17|| seems to be accidental. 
Here, we have shown that the susceptibility data for in- 
commensurate Peierls chains can be explained by a non- 
perturbative treatment of classical phase fluctuations. 
The exponential variation with temperature is due to the 
fact that in the pseudogap regime — ln% cx T _1 . Keeping 
in mind that our model is strictly one-dimensional and 
ignores amplitude fluctuations (which become important 
at temperatures of order T^), our theoretical curve for 
X(T) shown in Fig. agrees reasonably well with the 
susceptibility data ]17| . 

In summary, we have presented exact results for the 
average DOS, the mean localization length, the suscepti- 
bility and the low-temperature thermodynamics of dis- 
ordered incommensurate Peierls chains in the pseudo- 
gap regime, where only phase fluctuations are impor- 
tant. In particular, we have derived the exact frequency- 
dependence of p(tjj) which can be measured by means of 
angular integrated photoemission; we predict that at low 
temperatures p(u>) exhibits a maximum at lu = A s , the 
hight of which scales as T -1 / 3 . 
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FIG. 1. Frequency-dependence of the DOS given in Eq. 
(12) for g = 4A s £/v F = 0.4, 1.2,4.0, 12,40 and oo. 
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FIG. 2. The solid line is a graph of the DOS p(0) at the 
Fermi energy for classical phase fluctuations as a function of 
1/g = «f/4A s £. For a comparison, the dashed line shows the 
result found in the leading order Born approximation [5] and 
the diamonds give the DOS for Gaussian statistics [9,13]. 
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FIG. 3. Graph of the inverse localization length for 
g = 4A s £/v F = 0.4, 1.2, 4.0, 12, 40 and oo. 
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FIG. 4. Graph of the susceptibility x(T) calculated for 
£(T) = n s (T)/2m*T with n 3 (T) given in Ref. [16] and A S (T) 
determined by minimizing a generalized Ginzburg-Landau 
functional. The symbols represent susceptibility data from 
Ref. [17]. The dashed line is the DOS at the Fermi energy. 
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